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Abstract 




The Fokker-Planck coefficient for pitch angle scattering, 
appropriate for cosmic rays in homogeneous, stationary magnetic turbu- 
lence, is computed from first principles . No assumptions are made 
concerning any special statistical symmetries the random field may have. 
This result can be used to compute the parallel diffusion coefficient 
for high energy cosmic rays moving in strong turbulence, or low energy 
cosmic rays moving in weak turbulence. Because of the generality of 
the magnetic turbulence which is allowed in this calculation, special 
interplanetary magnetic field features such as discontinuities, or 
particular wave modes, can be included rigorously. The reduction of 
this result to previously available expressions for the pitch angle 
scattering coefficient in random field models with special symmetries 
is discussed. 

The general existance of a Dirac delta function in the pitch angle 
scattering coefficient is demonstrated. It is proved in this paper 
that this delta function is the Fokker-Planck prediction for pitch 
angle scattering due to mirroring in the magnetic field. The conditions 
under which this delta function contributes to pitch angle scattering 
are determined, and shown to be identical to the conditions under which 
first order mirroring occurs in the random field. These conditions 
are generally fulfilled in interplanetary and probably interstellar 
space. The implications of the delta function for the validity of the 
Fokker-Planck approximation are discussed. 
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I. Introduction 

The parallel diffusion coefficient for low energy cosmic rays in 
a random magnetic field has often been calculated from the small gyro- 
radius (or guiding center) approximation to the Fokker-Planck pitch 
angle scattering coefficient (Jokipii, 1966, 1967, 1968, 1971; Hall 
and Sturrock, 1967; Hasselmann and Wibberenz, 1968). But Klimas and 
Sandri (1973b), in a numerical calculation using the exact pitch angle 
scattering coefficient, have computed a parallel diffusion coefficient 
which differed markedly from that predicted using the small gyro-radius 
approximation. This calculation was limited to the special 
case of statistically isotropic magnetic turbulence with a Gaussian 
correlation function (and thus, a Gaussian power spectrum). Because 
of this limitation, the significance of this discrepancy was not fully 
appreciated because 1) in any case the small gyro-radius approximation 
was not expected to be accurate for Gaussian power spectra, (Jokipii, 
1971) but was still proposed as an accurate approximation for the more 
familiar power-law, power spectra, and 2) the source of the discrepancy, 
and especially its physical interpretation, was not dxscernable in the 
numerical computation. 

Recently, Fisk et al. (1974), and Klimas and Sandri (1973c), still 
within the framework of isotropic turbulence, but for power- law 
power spectra, calculated the pitch angle scattering coefficient 
without making the small gyro-radius approximation. They found that 
the exact result differed significantly from the approximate one, 
especially for pitch angles, 0 (measured relative to the mean field) 
near 90° where the Fokker-Planck coefficient contains a Dirac delta- 
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function in p = cos 6* The discrepancy found by Klimas and Sandri 
was then completely explained. For low energy cosmic rays in the 
Gaussian power spectrum, the delta function dominated all other 
contributions to the parallel diffusion coefficient. Furthermore, 

Klimas and Sandri (1973c) in calculating the parallel diffusion 
coefficient, showed that for power-law spectra, with any reasonable 
spectral index, the contribution of the delta function still dominated 
the contribution of the small gyro-radius approximation. For these 
spectra, Fisk et al. (1974) through a numerical calculation found that for 
p ^ 0, the pitch angle scattering coefficient is overestimated by the 
small gyro-radius approximation. Thus, the delta function is essentially 
the sole contributor to the Fokker-Planck estimate of the parallel 
diffusion coefficient in realistic models of the interplanetary magnetic 
field. A physical interpretation of the delta function, as well as a 
re-examination of the validity of the entire Fokker-Planck procedure 
in light of its existence, becomes necessary. 

Concurrent with the above developments, a number of papers (Jones, 
Birmingham, and Kaiser, 1973; Jones, Kaiser, and Birmingham, 1973; 

Kaiser, 1973; Jones and Birmingham, 1974; Volk, 1973; Volk, Morfill, 

Alpers and Lee, 1974) were presented in which it was argued that the 
Fokker-Planck formalism fails to correctly describe pitch angle 
scattering near 9 = 90° (p = 0). For mathematical simplicity, these papers 
have invariably discussed an idealized model of the random magnetic 
field in which the 6(p-) does not appear. We will demonstrate in this 
paper that the delta function is the rule, rather than the exception, 
and, therefore, must be taken 
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into account in studying the apparent failure of the Fokker-Planck 
formalism near p, = 0. 

In this paper we derive the Fokker-Planck pitch angle scattering 
coefficient with no restrictions on the statistical behavior of the 
random magnetic field, other than the assumption that it is homogeneous 
and stationary over several correlation lengths. Our calculation is 
baaed on the Liouville equation for the cosmic ray distribution function. 

We apply the quasi-linear and adiabatic approximations which, together, give 
us a Fokker-Planck equation for the mean cosmic ray distribution function. 
The mean distribution function is assumed gyrotropic in order to focus • 
our attention on pitch angle scattering. This derivation is done in a 
way which clearly demonstrates the existence of the 6(p) contribution to 
pitch angle scattering. But, the generality of our derivation enables 
us to study the conditions under which the strength of the delta function 
is zero or non-zero. In this way, we are able to determine that the delta 
function part of the pitch angle scattering coefficient is the Fokker- 
Planck prediction for the contribution of mirroring (Alfven and 
Falthammar, 1963; Northrop, 1963) to pitch angle scattering. It has been 
thought that the effects of mirroring were not contained, a priori, in the 
Fokker-Planck formalism, but had to be developed separately. Work along 
these lines has been done by Noerdlinger (1968), Quenby, et al. (1970), 
Cesarsky and Kulsrud (1972), and Jokipii (1973). In this last reference, 
Jokipii concludes that mirroring is not a significant pitch-angle scattering 
mechanism for an isotropic particle distribution. We show below, at least 
for isotropic magnetic turbulence, that the effects of mirroring dominate 
the spatial diffusion of low energy cosmic rays within the Fokker-Planck 


formalism. 
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In general, the interplanetary or interstellar magnetic field will 
have fluctuations in field magnitude which will produce particle 
mirroring to first order in the weak random field strength. For the 
isotropic models of the random magnetic field discussed above, the 
Fokker-Planck formalism predicts that mirroring is the most important 
pitch angle scattering mechanism for the parallel diffusion of cosmic 
rays. On the other hand, the existence of the delta function makes the 
Fokker-Planck formalism suspect. We feel that the isotropic field model 
calculation indicates that mirroring cannot be neglected, nor even 
included in a casual manner. In particular, further efforts at constructing 
a correct description of the behavior of those cosmic ray particles with 
pitch angles near 90° should probably not be based on idealized models of 
the random magnetic field which do not allow first order mirroring. 

II. Pitch- Angle Scattering Coefficient: Arbitrary Homogeneous Turbulence 

We have been able to construct the Fokker-Planck pitch angle scattering 
coefficient using two independent methods. The first method, which we 
will present here,, follows directly from first principles; i.e., the 
Liouville equation for the cosmic ray particle distribution function. 

The second method, outlined in Appendix A, is based on the standard 
Fokker-Planck approach (Chandrasekhar, 1943; Jokipii, 1971; Hasselman 
and Wibberenz, 1968). Although these methods are similar in some respects, 
they also contain surprisingly disimilar reasoning. The results of both 
methods are identical. 

Klimas and Sandri (1973a) have shown that, starting from the Liouville 
equation for the cosmic ray distribution function, the following truncated 
master equation for the mean cosmic ray distribution function can be 



constructed: 


+ oJ£f + eZ£ = (e 1 ) 2 <£'G £’> f (1) 

The right hand side of this equation represents the leading significant 
term in an infinite series expansion in the small parameter, e'. By, 
truncating this expansion at this point, we, in effect, make the familiar 
quasi-linear approximation. As a consequence of this truncation, we 
approximate the actual trajectory of a particle in the magnetic field, 
by the helical traiectory it would have in the mean magnetic field. We 
assume this helical trajectory for how ever long it takes the particle to 
travel approximately one correlation length along the mean field. Particles 
with pitch angles near 90° may take a long time to move this distance, 
however, and in this time the assumed helical approximation to the actual 
particle motion becomes suspect. It is this point which has lead to recent 
modifications of the quasi-linear approximation in the region of particle 
phase-space near, p, = 0. Our point of view, here, is to determine the 
actual predictions of the quasi-linear theory, which we may then use, in 
confidence, to further investigate the apparent failure of this theory 
and its possible modifications. 

Equation (1) has been written in a dimensionless form. The time 
variable is given by, t = t (v/X ), where t is the dimensional time, v is 
the particle speed, and is the correlation length in the random field. 

The term aiKf, is non-zero only if the distribution function depends on 
spatial position. Our purpose, here, is to fully investigate the pitch 
angle scattering coefficient that this theory predicts. Consequently, we 
can drop the term, oiKf, by assuming f independent of position. 
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The parameters, e and e* are defined by, 


e = 


\ <B> B’ 

-v 


( 2 ) 


where <B> and B' are, the mean and the rms magnetic field strengths 
respectively, and where P is the particle rigidity. We note that Klimas 
and Sandri (1973a) have given a rescaling of equation (1) which is necessary 
when e ^ 1. The version given here is appropriate when e £ 1. However, 
the results of the calculation we will present here can be applied to any 
range in e, so long as we remember that we must nave e' « 1 when e £ 1, 
but when e s 1, we must have T| « 1 instead . 

The differential operators, £ and £', generate the effects on the 
distribution function of the Lorentz force on a charged particle in the 
mean and random components of the magnetic ■ field, respectively. 

Convenient representations of £ and £’ , can be given in terms of the spherical 
coordinate system variables, 0, the polar, or pitch angle between the 
particle momentum and the direction of the mean magnetic field, and, 0, 
the azimuth, or phase ..angle of the momentum vector. With the definition 
of p, = cos 0, we are able to express these operators as, 

B 


£ = - 


B0 


and, 


s' - ^ (E-a-fi*) avV 1 


(3) 


(4a) 


or. 


s' = - ©'•«•£) ^ + aW 1 (&'•£•£) 

where, £ is a unit vector in the momentum direction, ^ is a unit vector 

in the direction of the mean field, 8’ = .B , /B' , and the tensors, 

— rms ' 
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Q and n are defined by, 


n ij = e ijk 


and, 


p.p. 

il 

n. . = 5 . . 7T- 

1J iJ 2 


(5) 


( 6 ) 


A full description of the integral operatior, G^, has been given by 
Klimas and Sandri (1973a). This operator, when operating on an arbitrary 
function of position, momentum, and time, operates as follows: 


ga Q 5’£> t ) = X ( e »^) A 


(7a) 


= JdX A (x(X), £ (X), t-X) (7b) 

o 

The streaming operator, J (e,X) , in equation (7a), shifts the phase space 

<6 

position (x,p) to (x(X) , j>(X) ) in equation (7b) which is along the helical 
particle trajectory in the mean magnetic field with (x, j>) for the starting 
point; i.e., 

x(\) = x - _r(X) (8) 

where 

W* A 

(9) 


r(X) = ^(e,X) 


A 

£ 


and 


A 

£ 


(10) 


(ID 


£(X) =C*(,?oM 
where 

C (e,X) = P + N cos eX - Q sin eX 
and, 

£ (e,X) = PX + ~ [n sin eX + Q(cos eX-1)] (12) 

The skew- symmetric tensor, Q, is defined above in (5) and the parallel 
and normal projection operators, P and N, are defined through, 
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and 


P. . = P. 

ij 1 K j 


(13) 


N.. = 6., - P.. 

ij ij iJ 

Now, we can rewrite equation (1) as. 


(14) 


i - * § ■ < e ’> 2< C i ^ «*»•& (H* 2 )’ 1 ] • 

(15) 

•J' T <wj, # (e,\) [-(£’• 0.£) (lVVh&’-n-g) > f(M,T-X) 

The integro-differential operator on the right side of equation (15) is 
quite complicated, as it acts on an arbitrary f(|i,0 s T-X). For the 
purpose of computing pitch angle scattering only, we assume a gyrotropic 
distribution function; i.e., we assume f independent of 0. Then, after 
averaging equation (15) over 0, we obtain, 


5f 

St 


(s') 2 !f D (t)*M 
2 S|i M- dp. 


(16) 


where, 

1 ,. T a A a 

V T> = ■ ~ J dX fa <(£-S*g. , (x)i (£' QO .0-£)> (17) 

p o o 

and where the star (*) in equation (16) indicates the convolution 

integral operation on the function of time to the right. Notice that 

this term, because it is assumed to be a function of p, and t only, is 

not affected by the streaming operator in D (t) . 

M 1 

We proceed by introducing the Fourier integral transform representation 
°f then allowing the streaming operator to operate on the 

resulting explicit dependence on the phase space variables. We obtain. 
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, ' 1 f & , 1 s3 r* -.3, r ,3. , i k.x 

V T) = " n J dX J d0 2tt ) j d k J d k e 


e 1 ~ * “ (x)) cl-a- < g/ (k) £/ &•)> •£• (M -£ 


( 18 ) 


W£th the assumption of homogeneous magnetic turbulence we can introduce, 

< £' (k) £' <£)> = (2n) 3/2 R(k) 6(k + k') (!9) 

where R (k) is the Fourier integral transform of the correlation tensor, 

R(r) - < £' (x) g/ (x + r) > (2°) 

No other assumptions concerning R (k) , or R (r) , will be made. 

Through the assumption of homogeneity we are able to rewrite the 
exponential terms appearing in equation (18) as, 

i k-x i k’ . (x-r(\)) -i k*r(X) /0 , x 

where £ (\) is given by equations (9) and (12). This result can be further 
rewritten to bring out the explicit dependence on the phase angle, 0. 

Then, the averaging over phase can be done. We will proceed in a particular 
coordinate system. Since D (t) is a scalar quantity, our choice of 
coordinate system is irrelevant. The end result of this calculation will, 
however, be written in an invariant form which then will apply in any 
coordinate system. 

In a cartesian coordinate system denoted by (1,2,3), we let the mean 
field be in the 3-direction. Then, both £ and k have parallel and 

A A A A 

perpendicular components denoted, for example, by (£)^=P»£ and (£) =N*£ 

with similar definitions for k and k . The magnitudes of these vector 

A a 2 1/2 

components will be denoted by |(p) | = p, and ] (g) \ - p 1 = (1-p, ) , and 

by k and k . The cartesian components of these vectors are defined through, 
tl 1 
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P x = u* COS 0 
p = V jl' sin 0 


and 


k_ = 


= 


k cos 
l 


k sin ^ 
1 


(22) 


P 3 = u k 3 = k 

Now, we can rewrite equation (21) as follows: 


-i k»r(X) -i k .r (X) -i k *r (X) 
e e e - ± 


(23) 


-i k nX » 
e \y S e 

m, n=-<» 


i(m-n) (^-0) e -iemX ( kj - U ‘ \ 

' m v ® ; n v e ' 


where the J's are the Bessel functions of the first kind. On substituting 
equations (19) and (23) into equation (18), we find. 


V T > = - i r i " 3/2 J dX I d3k e_i k " uX iji ( r n)(t ' 0) 


-i emX _ 


o o 

i 


(24) 


J_(— > J (-^-) t£-£?-J5Cfe) rfl- GtW -Jii 

m _ n e — - - ^ ■*“ 


Then, by introducing the tensor elements of R(k) , denoted by R-^Ck) , 

R 22 (~) > etc « » we are able to exhibit explicitly the phase angle dependence 

of the square bracket in equation 24. We find 
H- 2 *i 

[ 1 = — “ [-Q^-R^cos eX + (R^ R 22 )cos(eX+20) + 


+ (^12 +R 21) sin (eX+20) - (R^-R^) sin eX] 
From our earlier assumption of homogenity, we find 

RijCfe) = 

and, since the magnetic field is real 

R . .(-k) = R*.(k) 

ij ij 


(25) 


(26) 


(27) 


where the star (*) denotes complex conjugation. We will use these 
symmetry properties of the tensor elements in the following development. 
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Upon introducing equation (25) into equation (24), we explicitly 
exhibit all 0 - dependence in the integrand of equation (24), The 
phase angle averaging can now be done. Because of the orthogonality of 
the trigonometric functions in the interval, 0-2tt, only certain values 
of n, relative to m, in the infinite summation will contribute, thereby 
reducing the double infinite sum to a single one. The result is, 

WO " (^>' 3/2 (i ^ ) M JA • 

o 

• J m (~ ') [(R L 1 +R 22 )(e leX +e' :Le S-i(R 12 -R 21 Ke 15 :X -e" le?l )] + 


- el(2 ^ )+ e ' i(2i+eX)] + 

At this point, we must remember that D (t) is actually an integral 

operator (see equations (16) and (17). Further simplifications of this 

expression for D (t) depend on the function on which D (t) operates. We 
Li H 

must make the adiabatic approximation at this point to obtain the Fokker- 
Planck pitch angle scattering coefficient, even though there are good 
reasons for suspecting that this approximation may not be valid (Klimas 
and Sandri, 1971, 1973a). Klimas and Sandri (1973b) have shown that in 
the special case of isotropic, homogeneous magnetic turbulence, with a 
cosmic ray distribution finction which is essentially isotropic, the 
adiabatic approximation to the parallel transport of the cosmic ray 
particles can be formally constructed, but its accuracy remains in doubt. 
In the general case being considered here, no formal justification of 
the existence of the adiabatic approximation is available. However, 
because it is our intention to compute the Fokker-Planck coefficient as 
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a basis for further arguments on its validity, we will proceed anyway* 
The adiabatic approximation to equation (16) is the Fokker -Planck 
equation for pitch-angle diffusion, which can be written as, 


M_ll , <(Au) 2 > , M 

dT 2 du k At ; 


(29) 


where, 


= (e') 2 D (»> 
U 


At 


(30) 


In Appendix B we show that, by using the symmetry properties of equations 
(26) through (28), we find. 


«A i) g > = (€') g (i-u 2 ) r< j3 k « 
At 2/2tt J «=- 


J (^) 
» 6 ( em+k p,) m-l v e 


•f( R ll +R 22 ) J m-1 (^) - < R lf R 22> C ° S 21 + 

^ R 12 +R 21 ) J m+1 ( ~~e — sin + 1 ( r 12 “ r 21 ) •J m , 1 (“—)} 


(31) 


This result is applicable in a coordinate system with the 3-direction in 

the direction of the mean field. However, by using the relationships 

between the angle, \J/, and the components of k exhibited inequation (22), 

as well as the recursion relationships, J , (z) + J ,(z) = (— ) J (z) 

r * m-l N ' m+1' ' z ' m N ' 

and J . (z) - J . ( z ) - 2j'(z), we are able to rewrite equation (31) as, 
m— i m+i m 

.x2 


<i^T> m r d 3 k • « (em . k ^ . 

At /2=T J k ^ (em+k ll) 


k.u , 1 k.jjb’ k,,.R(k).k|. 

f-J 1 1 > T ( 1 ) r ^ Il. c.-H, 1 + 

1 nri-1^ e J m-l { e J L V 2 J 


+ <fV>" V Tr ^ ^ + 


2 „ 2 kiM*' 


k^V . .k,u, ' 

m r "i ' m ’ e 


-i(fV) J„ £*=-) J1 ££-■) Tr [C2.R (k)]} 

K j_M, 
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which now may be applied in any coordinate system. (The symbols, 

**Tr []**, indicate that the trace of the tensor inside the square bracket 
should be taken) . Equation (32) can be viewed as the pitch angle 
scattering coefficient which results from resonant interactions between 
the particles and waves in the magnetic field with wave numbers 
kjj = — . Similar expressions of this scattering coefficient have 
appeared previously for special forms of R (k) (Hasselman and Wibberenz, 
1968; Volk, 1973; Fisk, et al., (1974). In the next section we will 
consider the reduction of this general expression of the pitch angle 
scattering coefficient for any homogeneous magnetic turbulence to the 
previously available special cases. 

Ill . Isotropic, and Slab, Random Field Models 
i) Isotropic model for magnetic turbulence 

The magnetic turbulence is statistically isotropic if the tensor, 

R (k) , has the form, 

R(k) = R(k) (I - ) (33) 

where I is the unit matrix, and R (k) is a arbitrary scalar function of 
k = (k| which can, however, be related to the 1-dimensional power 
spectrum for the field components parallel to the displacement vector 
(Klimas and Sandri, 1973a) through the relation, 

R < k > - 1 k -k w 


By substituting equation (34) into equation (32), we find, for isotropic 
turbulence, 

<i^> _ j*d 3 k i = . fe <ly +e m) R(k) • (35) 
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This result, written in a dimensional form, has been given by Fisk, et al . 
(1974). 


All terms in equation (35) which correspond to m ^ 0 do not 
contribute to <(Ap,) 2 >/AT at \j, = 0. Notice that no part of the integrand 
depends on Thus , the integrations over \| / and k |( can be carried out. 
Then, 


<(Au) £ > 

At 



mtf>= K k x *(Vk 2 + (f!) 2 ) 


x M. 


^ + (^r) 2 4 c¥)0*(^) a )J t '_ 


(36) 

k' ' J k„= 


For a continuously differentiable magnetic field which is statistically 
homogeneous, R (k) must approach zero, as k approaches infinity, at 
least as fast as k 4 (Erdelyi, 1956; Batchelov, I960). Therefore, from 
equation (36),. <(Ah) 2 >/At = 0 at p, = 0. 

The m = 0 term in equation (35) contributes to <(A^) 2 >/At at p, = 0 
only. In fact, 



em 

M* 


Thus, it is the m ~ 0 term which contributes the 6(|i) part of the pitch 
angle scattering coefficient. Looking back at equation (32), we can 
see that this conclusion is completely general. 

The infinite series of resonant terms (m * 0) has been investigated 
by Fisk, et al. (1974), for power-law, pcwer -spectra. They found that, 
for Is E s 30, the entire series can be very well approximated by p, times 
the small gyro- radius approximation given by Jokipii (1971). This 
result applies over a wide range of values of the power law spectral 


index. 
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Klimas and Sandri (1973b, 1973c) have shown that the m = 0 terra 
dominates the parallel spatial diffusion coefficient for both power 
law and gaussian power spectra for low energy particles (e»l). The 
parallel diffusion coefficient is given by Klimas and Sandri (1973b) 


as 




(38) 


for T)«| and e»|, where Xp is the parallel integral length (essentially 
the correlation length in the random field). The quantity, A ^ (e ) , is 
related to the pitch angle scattering coefficient through (Fisk, et al. 
1974), 




(39) 


For e»l, equation (37) can be simplified considerably, Using 
^1 (-— ) "et ^ (~) , and equation (34), we find 


CO 

<(Au) g > | _ 6<H)T| fdk kP (k) 

At m=0 etc ^ *t) 


(40) 


On substituting equation (40) into equations (38) and (39) , we reproduce 
the earlier results of Klimas and Sandri (1973b, 1973c) as well as 
obtain an excellent approximation to the parallel diffusion coefficient, 
as predicted by Fokker-Planck theory, for low energy cosmic rays. A 

general feature of this result (in fact, for the S(^) contribution to 

2 

<(Ap) >/At in any kind of homogeneous turbulence) is that the term 
2 -1 

(e A.,(e)) , in equation (38), is independent of e, , for e»l. Thus 


K„ can be written, 


K li = 3 v Sbp 


where the mean free path, which is given by 


(41) 


is independent of regidity, 


( 42 ) 
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ii) Slab model for magnetic turbulence 
In the slab model we let 

R(k)= N R(k t] ) 8(k x ) (43) 

where R(kj|) is an arbitrary even function of k^. Thus, the random 
field components are normal to the mean field, and they are functions 
only of distance along the mean field. The correlations in the random 
field are assumed cylindrically symmetric about the mean field. 

This model is closely related to others which have been studied 
previously. In Jokipii's (1972) plane polarized model, the random field 
is assumed to have only one orthogonal component, We will see, shortly, 
that the pitch angle scattering coefficient computed here is simply twice 
Jokipii's. Jones, Birmingham and Kaiser (1973) have used the plane 
polarized field model to investigate modifications of quasi-linear theory 
near p, = 0. Volk (1974) has also studied such modifications in a random 
field made up of transverse Alfven waves. A common feature of all of these 
models is that the 6(p.) does not appear in the Fokker-Planck pitch angle 
scattering coefficient. 

With the choice for R(k) given in equation (43), we find, from 


equation (32), that, 


. 


„ g dlgbU . J d \ 6(V R(k||) ack^+em) 


<„■=, ) 2 j 2 Aal) 


(44) 




m 


We see immediately that there can be no 5(p.) contribution to 

2 

<(Aji) >/At since the m = 0 term in equation 44 is zero. In fact, because 
of the 6(k ± ), only the m = +1 terms are non-zero. All other terms in the 

ero k ^ 1 

expansion are zero because ( ■ -) J (-^ — ) - 0 as k - 0, and the integrand 

in 
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can only contribute at = 0. Thus 

“ (• , ) a ^Tt <i?t> R( fc ,i = p <«> 

(Note: 6(k ± ) is normalized to 2rr instead of 1). This result is exactly 
twice that of Jokipii's (1972) plane polarized field model calculation 
in which there is only one orthogonal field component instead of two. 

On substituting equation (45) into equation (39), we find, 

*. ■ '{1 £ * «-* 2 > *<;> 

■ 3 Vi J’ ¥ <‘-<!> 2 > «■> («> 

e 

which is identical to the small gyro-radius approximation to A,| that 
Klimas and Sandri (1973a) have computed in their isotropic model. This 
correspondence between the small gyro-radius approximation in isotropic 
turbulence and the exact result in the slab model was first pointed 
out by Jokipii (1971) . 

IV. Mirroring and the Delta Function 

The fact that the 6(p,) contribution to <(Am.) 2 >/At vanished in the 
slab model was one of the first strong indications that this contribution 
represents mirroring in the random magnetic field. The guiding center 
approximation (Alfven and Falthammar, 1963; Northrop, 1963) to the motion 
of charged particles in the random magnetic field applies to particles 
whose motion along the field is so slow that they can resonate only with 
very short wavelengths where the power density is typically negligible. 
This could be the situation for a very low energy particle, or a higher 
energy particle moving in a weakly turbulent field with pitch angle very 
near 90°. The latter situation is the one being studied here. 
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Within the guiding center approximation, the pitch angle changes 
with time according to, 

- M e • [(e-V) B] (47) 

12 A 

where, M = — /Bm, is the magnetic moment, and e is a unit vector in 
the direction of the local magnetic field. Of course, the pitch angle 
that enters into this expression is the angle between the direction of 
the particle momentum, and the local exact field (actually the angle 
averaged over a gyro-period). In Appendix C we demonstrate that, to 
first order in 1 ] in a weakly turbulent magnetic field, this pitch angle 
is identical to the pitch angle relative- to the mean field averaged over 
a gyro-period By introducing B = <B> + B-J we can rewrite equation(47) . 
To first order in T\, we find, 

^ \ a-\£) T) J3 • [<£ . V) £■ (k>] (48) 

where 1 / 2 ( 1 -p )T) is the magnetic moment, which we consider a constant 
of the motion in this guiding center approximation. Since the pitch 

angle scattering coefficient is quadratic in Ap, we see that, in order 

2 

to obtain an Q (p ) contribution to that coefficient, we must find an 

0 (T|) contribution to d^/dr through equation (48). Equation (48) is 

non-zero only if the random field contains a component in the direction 

of the mean field which also has spatial gradients in that direction. 

In the slab model, (3 .|3 1 (x) = 0, and there is also no 6 (p) contribution 
2 

to <(Ap) >/At. In this section we will show that if, and only if , 
equation (48) is non- zero, do we find a 5 (p.) in the pitch angle 
scattering coefficient. Thus, only when there are 0 (11) changes in p 
due to guiding center motion, is there a 5 (y.) in <(Ap,)^>/Ar. In view 
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p 

of the fact that the 5 (jj,) is the only part of <(Ap) >/At which is 
non- zero at p = 0, and therefore, is the only part which can possibly 
reverse the parallel motion of the particles, we ascribe the 6 (p) to 
mirroring. The fact that mirroring exhibits itself as a 5 (p) in 
<W >/At can be understood from the asymptotic nature of the Fokker- 
Planck equation. This equation is asymptotic in small T], and, in an 
infinitesimal random field, only particles very near p, = 0 could possibly 
be mirrored. In fact, the delta function is probably a crude approximation 
to a sharply peaked function of p. 

It is interesting to note the results of a calculation which we 

p 

present in Appendix D. There, we compute <(Ap) >/At from equation (48) 
using the Fokker-Planck formalism. For the position of the particle, 
x, along its trajectory in the field, as it appears in p'(x), we substitute 
the position of its guiding center. This approximation to the position 
of the particle is exact in the limit of zero particle energy (e = ») . 

The result of this calculation is identical to the result which we will 
present in this section when we study the m = 0 term (the 6 (p) contribution) 
of equation (32) in the limit, e = «. Thus, the Fokker-Planck coefficient 
for pitch angle scattering, computed in the guiding center approximation 
to the particle trajectory, also contains a 6 (p) term which is in exact 
agreement with the 6 (p.) term obtained from the quasi-1 inear approximation 
to the particle trajectory, in the limit, 0 = 0 =. 

From equation (32), with m = 0, we find, 
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A necessary property, of R (k) is that it be non-negative i.e., the 
quadratic form, 

V * R • V*s 0 (50) 

for any complex vector, V (Batchelor, 1960). In particular, k * R*k ( | 

must be either positive or zero for any k. Notice also, that the rest 

of the quantities in the integrand of equation (49) are non-negative. 

Therefore, barring special situations which we will discuss in a moment, 

R (k) , must be zero everywhere in k in order for the 6 (p) to not 

contribute to pitch angle scattering. But then, R (r) is zero for 

II > II 

all from which we can conclude that there is no random field component 

parallel to the mean field (slab or plane polarized field model). With 

no random field parallel to , we see from equation (48) that mirroring 

is impossible. This argument can be carried in the reverse order. It 

is clear that if equation (48) indicates no mirroring because (3*(3 1 (x) 

is zero for all x, then R (r) = 0, and finally R (k) = 0 for all k. 

1MI ~ IM! ~ 

2 

In this case there is no 6 (p) contribution to <(Ap) >/At. 

There is a special situation in which R^ ^ (k) ^ 0 everywhere in k, 
and yet the integral in equation (49) is zero, and <(Ap)^>/AT does not 

contain a 6 (p) . It is possible that R (k) is non-zero only at point 

II > II 

in k where other terms in the integrand are zero. The Bessel -function 
cotains isolated zero’s which, nevertheless we will not consider because 
the positions of these zeros are rigidity dependent through e, and R (k) 
does not depend on rigidity. On the other hand, we could have R (k) 
non-zero only when k^ =0 without contributing to the integral. But, 
we still have, 

k il R ||,|l® = °- 


(51) 
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This statement is the Fourier transform of 




= 0 


( 52 ) 


Barring non-dif f erentiable functions of the type discussed by Wiener 
(1933), we can conclude that the random field component in the direction 


of the mean field is independent of position along the mean field; i.e., 
(P-V)(|.- J3 '(x» = 4t4-v)p'(x)] = 0. But this is exactly the condition 
that equation (48) be zero. This argument can also be reversed. Starting 
from a zero for equation (48), we can conclude equation (52), then 
equation (51), and then we can conclude no 6 (p) in <(Ap) 2 >/At. In general, 
a random field component along the mean field is allowed, but this component 
must be a constant along the field; in the slab model the constant is zero. 

Fr om the above arguments, we conclude that the 5 (il) in <(AlQ 2 >/At 
represents the physical phenomenon of mirroring in the random magnetic 
field. 


In isotropic turbulence, we have seen, when e » 1 (low energy 

O 

particles) the 6 (p) contribution to <(Ap) >/At becomes independent of 
rigidity. From equation (49), we find in the same approximation. 


| ( (±-V?)\ 2 kitel r ri 3 k r k|| -R(k) .k|| 

AT m=o 2 /2rr J ] kii I 


(53) 


Thus, the 6 (p) part of <(Ap) 2 >/AT becomes rigidity independent, for 
low rigidities, in any homogeneous magnetic turbulence. For isotropic 
turbulence, Klimas and Sandri (1973b, c) have shown that this 6 (p) 


contribution dominates over all other scattering mechanisms. Thus the 
rigidity dependence of is determined by the effects of mirroring: 

SiFP '*’ S r ^-§£<lity independent. For arbitrary turbulence, however, 
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we do not know the relative strengths of the contributions of the 6(^) 
terms, and the non-5 (|i) terms, to the parallel diffusion coefficient. 
Thus, in contrast to the isotropic turbulence, we cannot determine the 
rigidity dependence of . 

As we mentioned above, equation (53) is identical to the result 
(given in Appendix D) of computing <(A|i,) 2 >/AT within the guiding center 
approximation to the particle motion. 

V. Conclusion 

We have constructed the Fokker-Planck coefficient for pitch angle 
scattering of cosmic rays in otherwise arbitrary, but statistically 
homogeneous, magnetic turbulence. Our result was obtained both from 
first principles, and through the Fokker-Planck formalism. The reduction 
of our expression to previously available scattering coefficients, 
calculated in special models of the random field, has been discussed. 

We have shown that the pitch angle scattering coefficient contains 
a Dirac delta function, 6 (p,) , in p, = cos 0, where 0 is the pitch angle. 
We have, further, proved that this delta function is the Fokker-Planck 
prediction for the contribution of mirroring to pitch angle scattering 
in a weakly turbulent magnetic field. 

The 6 (p) does not contribute to pitch angle scattering when, within 
the guiding center approximation to the particle motion, the pitch angle 
is a constant of the motion of the particle 0 (T)) in the random field 
strength. This condition is met when the vector component of the random 
field, which lies in the direction of the mean field, is independent of 
distance along the mean field. The slab, or plane polarized, or 
linearized Alfven wave models of the random field, all of which have 
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no random field component along the mean direction, fall within this 
class. Magnetosonic waves propagating across the mean field have a 
random parallel field component which is independent of distance along 
the mean field, and fall within this class. Typically, however, the 
observed interplanetary field does not fall within this class. We 
expect mirroring in the interplanetary field, and we expect a delta 
function in the appropriate Fokker-Planck pitch angle scattering 
coefficient . 

The fact that mirroring exhibits itslef as a 5 (p.) in the Fokker- 
Planck pitch angle scattering coefficient indicates that this particular 
pitch angle scattering mechanism is misordered within the Fokker-Planck 

formalism. In this formalism it is assumed that the effects of the 

2 

random field on the motion of the particles is 0(7] ) in the random 
field strength. When the effects, on the cosmic ray distribution 
function, or mirroring are computed correctly we will most certainly 
find them entering in a lower order than 0 (T] 2 ) , in a non-Fokker-Planck 
type equation. 

Other investigators have concerned themselves with the modifications 
of the Fokker-Planck formalism which are necessary even in the special 
random field models in which mirroring does not occur. Even though 
these modifications are extremely difficult to compute, they nevertheless, 
in principle, are minor modifications to the quasi-linear , or weak 
coupling approach which leads to the Fokker-Planck equation. The phenomenoi 
of mirroring has not been properly ordered within the approach, and the 
presence of the delta function in the pitch angle scattering coefficient 
presents, a serious obstacle to these modification schemes. In this paper 
we have demonstrated the existence of the delta function, as well as its 
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connection to mirroring, in order that further discussions of the 
validity and/or modifications of the Fokker-Planck approximation can 
be based on the actual Fokker-Planck predictions in realistic field 
models , 

We would like to acknowledge useful discussions with Dr. L, A. Fisk. 





Appendix A 

With our notation, the momentum equation of motion for a charged 
particle in a magnetic field, can be written as, 

A 

ffe = e [n + T| O'] . p (A. 1) 

dT 

where an ^ * € ijk^k* T h us > for p, = J3 • £ we find 

= € ' * £?'* 2) (A. 2) 

or, 

IS - e '<E • a • el') (a. 3> 

where, in obtaining equation (A. 2) , we have used = 0. By formally 
integrating equation (A. 3) we can obtain, 


<(Ap) 2 > = - (e 1 ) 2 JdX Jdx£(\).Q. <g/(X)£'(X’)> -O-pCX') 

o o (A. 4) 

where g/ (\)^3 ’ [x(X) ] , and x (X) and p (X) are the time dependent 

coordinates given by equations (8) thru (12) which describe the zeroth 

order helical approximation to the particle motion. Now, by introducing 

the Fourier integral transform representation of g* (x) , as well as the 

assumption of homogeneity, as expressed through equation (19), we find, 

2 At At 

<(A^) 2 > = JdX JdX Jd 3 k e ifc ll^ X - X V£±- ( ^ W -£L (X, >>. 

° u o 

(A. 5) 

We also choose a special coordinate system in which to proceed here. In 
a cartesian coordinate system, we let the mean field lie in the 3-direction, 


as in the main text, but we further rotate our coordinate system so that 
the perpendicular component of the momentum, (p ), lies In the 1-direction. 





Then, on introducing the expansion, given by equation (23), of the 
exponentials in equation (A. 5), we obtain. 


<(Ap,) 2 > = - -^Tp > ? dk l! e ik l^ (X_XI) rS l e i( m*-n) t|f 

(2tt) 3 / 2 J 0 J 0 i 11 J i 


ineX -imeX' 
e e 


J' dk X 


k.u,' 
j ) 
m e 


• [p (X).Q*R(k)*Q.p,(X')] (A.6 

By introducing the tensor components of R (k) , we find, 

0 / i \ 2 , - 2. At At » v 2rr . . . . 

<(A^) 2 > = J dX J dX > J dk e xfc t |P.(X-X ) j d+ e Km-n)t 


2(2rr) 


O O -oo 


ineX -imeX' 


fdlc, k, jn (i±2) J (i±tLl) . (A.7) 

_L x x a y m € 


.. C(R 1X +R 22 ) cos e(X-X') - (R^-R^) cos e(k+X') + 

- ( R i 2 +R 21^ sin e (X+X* ) - (R^-Rgj^) sin c (X-X 1 ) 1 


Now, following the usual Fokker-Planck arguments, we must find the parts 
of this expression which grow as At. Consider, for example, the integral, 

I x - Jd X Jdl' ) e in^ e -imeX' cog e(X . V) (A . 8) 

o o 


which is associated with the term (R^+R^) in equation (A.7). By 

introducing a new variable of integration S = X'-X, for X', and inverting 

the order of integration, we find, 

At At-S o At , 

I x = [Jds JdX + JdS Jd^ d’ l(k ll ,i+em)S cos e S e l(n - m)eX (A. 9) 
o o -At -S 

When n = m, we do obtain a term which grows linearly with At, for large 
At, given by. 


At tt [6(k |lM , + (nH-l)e) + 6(k ||M , + (m-l)e)] 


(A. 10) 
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A further examination of equation (A. 9) reveals no other terms with the 
linear growth property. In order to obtain the Fokker-Planck coefficient 
from equation (A. 7), we replace 1^, in equation (A. 7) with 

At n rr [6(kj,M. + (m+l)e) + 6 (k^ + (m-l)e)] (A. 11) 

This procedure can be continued for the rest of the time integrals, and 

the steps which follow are identical to the steps taken in the alternate 

2 

derivation of <(Ap.) >/At given in the main text. As we mentioned there, 
the results of the two calculations are identical. 

The point of this appendix is to note the difference between the 
arguments given here, and those given in the main text. From Equation (A. 9), 

we actually find a variety of terms, some of which do not grow unbounded 

2 

with At, some which grow linearly with At, and some which grow as (At) . 

From the other time integrals in equation (A. 7) we find similar results 
with, in addition, terms which grow as At times trigonometric functions 
of At; these terms oscillate with amplitudes which grow unbounded in time. 

In the face of these unbounded oscillations, the usual Fokker-Planck 

2 

argument for choosing the terms which grow as At in <(A|i) >, and dropping 
all others, becomes difficult to support. Even ignoring this problem, 
we still have the problem of choosing the At terms by arguing that At 
is large compared to the interaction time, and small compared to the 
relaxation time, when in: fact these two time scales are not clearly 
separated for particles in a magnetic field. 

In comparison, from equation (24), notice that we never face these 
problems in the derivation of <(A^l) 2 >/At given there. The inner time 
integration in equation (A. 7) is replaced by the phase averaging 
integration over 0 in equation (24). The effects of carrying- out the 
Fokker-Planck argument here, and the phase averaging there are identical; 



both procedures pick the same terms out of the infinite expansions for 
retention. Thus, the results of both approaches are identical, but, 
the reasoning contained in the Fokker-Planck approach is much more 
difficult to support. 
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Appendix B 

The adiabatic approximation to equation (28) permits the time 
integration to be done. The result is 


D («) = (2tt)" 3/2 Jd 3 k S J m • 

M* d 0 m=-® e 


(1^1) [ iCl^+Rgg) CC(e(l-m)-k|,ji) + £(-e(l+m) - k^p,)] + 


+ ( R i 2 ~^ 2 i) CC(e(l-ra)-k 1( p) - C(-e(l+m) “ k t |l*)3 } + 

-i(Rll- R 22) [J m . 2 (illil) e i2 * C(«(l-»)-k,l») + W— > 

e e 

-(Rl 2 + R 21 ) [J m . 2 (^) e 12 * C(«(l-»)-k,H) - 

e e 


(B.l) 

e" 12 * C(-e(l*Hn)-k,|^) 
e l2 ^ £(-e(l+m)-k,||i)] 


The function £(x) is the zeta function defined by (Heitler, 1954) 

T 

C(x) = - i lira J*dX e iXx = \ - in6(x) = - C*(-x) (B.2) 

T-*® O 

p 

where - is the principle value of 1/x. 

X 

In the terms that multiply £(e(l-m)-k||p,) , let m' = - m and k' = - k, 
and use the relation J_ m ( z ) = (-1)™ J m ( z ) a l° n g with the symmetry relations 
(26) and (27) and (A. 2) to rewrite (B.l) in the form 


2 

y-> = W~ 3/Z ^ I d3k i.» J m {V£-C*)U(Rii+R22> + 

( R 12- R 21>] " J tiri-2 (R ll' R 22 ) ^- e '*' C " e ' 

- J mf2 ( R l2 +R 2i^ ei2t £ + e_i2t £*13 (B.3) 

where J m = , C = £(e(l+m)+k ||M ,) , and £* = e*(e(l+m)+k ||M ,) . 

The terms containing zeta functions can be rewritten as 



COO - C*00 - - 2TTi 6 (x) 

e i2 ^£(x) " e i2 ^*(x) = 2i[sin 2^ (|) - tt cos 2f (6(x))] 
e l2 ^C(x) + e l2i C*( x ) = 2[cos 2ijr (|) - tt sin 2if (6 00)3 

(B. 4) 


Now consider terms of the form 


J d k 5 _ J m J 


m-co 


tiH-2 


( R 11 _R 22) Sin 2 ^ 

( R 12 +R 21) cos 21 f 


(B.5) 


e(l+m)+kjj^i 

and let n = m + 1, so that equation (B.5) becomes 


« 3 03 

\d k. 2 J -J 
J n=-a> n+1 


( R 11 _R 22) sin 2i|r 
^ R 12 +R 21) c0s 2 ^ 


en+k n 


(B.6) 


If we let n -* -n and k -* -1c in equation (B.6), then equation (B.6) is equal 
to minus itself and all terms in (B.3) which contain principle value 
contributions are identically equal to aero. Equation (31) follows 
immediately o 
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Just for the purpose of developing this argument, in this Appendix, 
we will introduce a special notation which is different from that con- 
tained in the rest of this paper. Let, 

M-W = • Z (C.l) 

be the cosine of the pitch angle relative to the local (exact) field. 
Within the guiding center approximation to the particle motion, we 
assume (B/B) constant over a gyroperiod, and assume that the particle 
moves in a' helical trajectory in the field. Thus, 

p.(T) = (|) • 5 + N cos st - Q sin ct] . £(o) (C.2) 

B 

where the projection matrices, P, N, and Q, are identical to those 
defined by equations (5), (13), and (14) except that they are based on 
the local, rather than the mean field. We now introduce the notation 
<h> = p . p, for the cosine of the pitch angle relative to the mean 
field. The quantity which enters equation (47) is \x (t), which is the 
average over one gyroperiod, of p, (t) . 

M- (t) = (p • P • £(o) s [x (C . 3) 

The time dependence of <p,> is given by, 

<p> = 4 * LjP + N cos er - 0. sin st] • £(o) (C.4) 

and its average over one gyroperiod is, 

<|j> = P .f + T l &' “I 

Li-l+2T|(£ «3 1 ) + Cg* -3* )1 1/2 J ^ 


(C.5) 
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To 0 (Tl) in the radom field strength, <p> and p are identical when . 
averaged over a gyroperiod. This averaging is implied in equation (47) 
since it follows from the guiding center approximation. In the rest 
of this paper we therefore, simply use the symbol \i to stand for the 
cosine of any of the relevant pitch angles. 
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By introducing the Fourier transform representation of the random 

field into equation (48), we find, 

3/2 


f; = - | d-M 2 )^) Jd 3 k |_.[(j|.V) e^^V (k)] 


(D.l) 


or, 


1 3/2 


? (1-/)^) Jd 3 k(|.p'(k))(|-k)e ~~ (T) ' 

2 

Now, on remembering that, in this notation, (1-p, )T)/2 is the magnetic 
moment of the particle which is an adiabatic invariant of the motion 
of the particle in the guiding center approximation, we can formally 
integrate equation (D.2) and form its square to find. 


<(A,) 2 > - - (i^) 2 (i) 3 Jd l Jdl-Jd\ jA- 


(D.3) 




If the magnetic turbulence is statistically homogeneous, we have 
equation (19), which can be inserted into equation (D.3), to obtain, 

«A,) 2 > - ( ( i ^) 2 <^) 3/£ Jdl Jdl' Jd\ e 1 

o o 

• fell • R(k) • Jk (| ] 

As e xplained in the main text, we approximate the position of the 
particle, by the position of its guiding center. Thus, k* (x(X)-x(X' ))= 
= s k||p,(X-X I ) . Following the usual Fokker-Planck argument (see Appendix B) 
for obtaining the part of equation (D.4) which grows linearly with At, 
when At is large, we obtain, 
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2_2 


<(AiO > _ r (1-yx 6(kL) r\,3, ,k„.R(k).k n 

At ” C 2 > ^ J d k 7, ^7- " ) 


l k l 


This result is identical to that given by equation (53) which was 
obtained, in the low energy limit, from the usual quasi-linear 
approximation to the particle motion. 


(0.5) 
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